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Abstract — For the last decade, it has been known that reflected 
GPS signals observed with specialized instruments could be used 
to measure sea level. In this letter, data from an existing geodetic- 
quality GPS site near Kachemak Bay, Alaska, are analyzed for 
a one-year time period. Daily sea-level variations are more than 
7 m. Tidal coefficients have been estimated and compared with 
coefficients estimated from records from a traditional tide gauge at 
Seldovia Harbor, ~30 km away. The GPS and Seldovia estimates 
of M 2 and S 2 coefficients agree to better than 2%; much of 
this residual can be attributed to true differences in the tide 
over 30 km as it propagates up Kachemak Bay. For daily mean 
sea levels the agreement is 2.3 cm. Because a standard geodetic 
GPS receiver/antenna is used, this GPS instrument can measure 
long-term sea-level changes in a stable terrestrial reference frame. 

Index Terms — GPS, multipath, reflections, tide gauge. 

I. Introduction 

T HE TECHNOLOGY for measuring sea level at a coastal 
site is well established. We are nearing almost two cen- 
turies of experience with self-recording tide gauges, with the 
most common one consisting of a float isolated from high- 
frequency wave motion by a stilling well [1], but acoustic, pres- 
sure, or microwave-based systems are also now in use [2], [3]. 
The conventional tide gauge measures relative sea levels, i.e., 
relative to benchmarks on the adjoining land, which is the 
quantity of interest for most applications such as navigation 
and coastal flooding. Other applications, particularly estimating 
changes in the global ocean volume and related aspects of 
the global sea-level budget, require absolute sea-level measure- 
ments, i.e., relative to the geocenter or to a fixed terrestrial refer- 
ence frame [4], [5]. Effects such as glacial isostatic adjustment, 
coseismic and postseismic deformation, and land subsidence 
make it challenging to use conventional tide gauges to measure 
absolute sea level. These kinds of “local” ground motion at tide- 
gauge stations can, however, be measured with continuously 
operating GPS receivers, and the resulting information can 
be used to convert between relative and absolute sea levels. 
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The critical need for such land-motion information has led to 
an international campaign to deploy GPS receivers (and other 
geodetic systems) at a wide array of globally distributed tide 
gauges [6], [7], 

It was recently demonstrated that geodetic-quality GPS in- 
struments installed to measure crustal deformation, but sited 
very near the coast, could be used to measure sea level directly 
and without local ground motion errors [8]. In this experimental 
setup, power measurements made by the receiver are used to 
deduce the distance between the GPS antenna and the sea, 
which is relative sea level, yet the GPS station also provides, 
by design, the data necessary for placing those estimated sea 
levels into a true terrestrial reference frame, i.e., into absolute 
sea level. This is distinct from previous GPS reflection studies 
that used either specialized instruments or installations (see, 
e.g., [9], [10], and references therein). 

In this letter, we examine a one-year data set from a chal- 
lenging GPS site in Kachemak Bay, Alaska. Tidal variations at 
this site are among the largest in the United States. We are able 
to estimate both high- and low-frequency sea-level variations, 
which can be compared with similar signals from a nearby 
conventional tide gauge. The large tidal variations present some 
new additional challenges for deducing sea level from GPS 
reflections, a topic discussed in detail in the next section. 
These data sets present a unique opportunity to investigate 
such measurements in a challenging, and clearly macrotidal, 
environment. 

II. GPS Data 

The GPS site at Peterson Bay (PBAY) is located on a small 
rocky island on the south side of Kachemak Bay, in the Cook 
Inlet of Alaska (see Fig. 1). Due to the dense vegetation cover 
of most of the region, only small rocky islands remain free of 
trees and thus have suitable sky view for GPS observations. 
PBAY was installed by the University of Alaska Fairbanks as 
part of a project to better understand crustal deformation in the 
area surrounding Kachemak Bay (see Fig. 1 and Appendix). 
The unit provides coverage of the south side of the bay, which 
has bedrock at the surface, to ensure that one could detect any 
contrast with the north side of the bay, which has sediments at 
the surface. 

Reflected GPS signals (multipath) are one of the major error 
sources in crustal deformation studies, where the objective is to 
measure the position of the GPS antenna phase center. Mul- 
tipath effects measured by geodetic-quality GPS instruments 
can also be used to accurately measure environmental signals 
(see [8] and references therein). A defining feature of these 
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Fig. 1. Photograph of PBAY GPS site at low and high tides. 

multipath studies is the use of signal-to-noise ratio (SNR) 
data instead of the carrier phase data used for positioning. 
SNR data are dominated by long-period variations due to the 
receiver-satellite distance and the gain pattern of the antenna. 
If these long-period variations at low elevation angles (< 30°) 
are removed with a low-order polynomial, one is left with only 
the multipath effect. 

For a horizontal planar surface, the model for GPS multipath 
is straightforward. The multipath effect is driven by the instan- 
taneous excess propagation delay d of a reflection with respect 
to the direct path 

d = 2Hsine. (1) 

H is the height of the antenna above the reflecting surface (or 
reflector height) and e is the elevation angle of the satellite with 
respect to the horizon [11]. When the surface is stationary, the 
instantaneous phase (j) = kd (where k = 2ir/\ is the wavenum- 
ber and A is the GPS carrier wavelength) is driven by a constant 
factor of AnH/X times the variable sine [11]. Since GPS 
multipath measures the distance between the antenna and the 
reflecting surface, the highest water levels produce the smallest 
reflector heights and the lowest modulation frequencies. Fig. 1 
shows representative SNR data from the PBAY site. The SNR 
data clearly show significantly different multipath frequencies; 
these are anticorrelated with sea-level heights measured by a 
tide gauge. 

If H is constant, spectral analysis of SNR data with sine 
as the independent argument yields a peak whose frequency 
is linearly related to H. However, this model is not accurate 
enough for the PBAY data where a new model appropriate for a 
moving planar surface is required. We start with the derivative 
of d with respect to time 

d = 2He cose + 2H sine. (2) 

Equation (2) can be interpreted as the sum of a dynamic satellite 
term and a dynamic surface term [12] . An angular delay rate can 



Fig. 2. GPS SNR data (converted from dB-Hz to volts/volts and detrended) 
for two GPS satellite tracks at high and low tides. 


Sea Level 
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Fig. 3. Time series of sea-level estimates from (circles) GPS and (line) tide- 
gauge data for Seldovia. 

then be defined 

° 7 dd dd / d sine d 

d = — = ^7 / *+ = ( 3 ) 

asine at / at ecose 

or 

°d= 2H + 2H tan e/e. (4) 

The first term on the right-hand side of (4) is the usual result 
for stationary surfaces discussed earlier, and the second term is 
a height-rate correction to the otherwise biased spectral peak 
height. 

We solve for the two unknowns iteratively. We first esti- 
mate H values with a Lomb-Scargle periodogram to obtain 
a preliminary height time series. We then estimate H and 
calculate tan e/e values to produce height corrections (see 
Fig. 2 and Appendix). These are removed from the preliminary 
time series, resulting in a final height time series (see Fig. 3). 

GPS data collected between June 26, 2011 and June 30, 
2012 were analyzed. The locations of the transmitting satellites 
were computed using the broadcast ephemeris; these are subse- 
quently used to compute the local azimuth and elevation angle 
(and its rate of change) of the satellite with respect to the GPS 
receiver. A total of 39 rising or setting GPS satellite tracks were 
used. Minimum and maximum elevation cutoffs were set to 5° 
and 25°. Above 25°, the multipath signature is small, and the 
reflecting region (Fresnel zone) often includes both land and 
water (see Fig. 3 and Appendix). The reflecting region varies 
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with elevation angle and reflector height, but in general, each 
track provides an average over thousands of square meters. 
For this particular site, an azimuth mask of 90° to 260° was 
imposed. This mask removed reflections that are impacted by 
the trees behind the GPS antenna. 

In order to report a reflector height, we required the use of 
800 SNR measurements (~14 min), a normalized spectrogram 
peak of 1 dB-Hz, and a peak-to-noise ratio greater than 3. “Non- 
physical” reflector height peaks were discarded, i.e., reflector 
heights greater than 1 1 m or less than 3 m. After these quality 
control measures were imposed, ^80 sea-level retrievals were 
produced per day. However, it must be emphasized that these 
retrievals are not evenly spaced in time as they would be for a 
dedicated tide gauge. 


III. Discussion 

The definitive test of the GPS -based sea level estimates 
is comparison with measurements made with standard tide 
gauges. As there is no exactly collocated tide gauge at the 
PBAY site, some allowance must be made for small, but real, 
ocean differences between sites. At the Seldovia tide-gauge site 
(30 km distant), a long multidecade time series of hourly sea- 
level data is available, including data concurrent with our GPS 
data. There are also a handful of other historical gauges in 
Kachemak Bay; their time series are all very short (29 d), but 
they are still useful for delineating general tidal characteristics 
within the bay. 

It is most useful to separate the comparisons into two 
regimes: low frequency and high frequency (once/day and 
higher). The high-frequency comparisons could be done with 
(for example) 6-min or hourly measurements, once the GPS 
estimates have been appropriately processed and resampled. 
The high-frequency regime, however, is so dominated by large 
tidal signals, and the corresponding sea-level differences are 
likely dominated by small changes in the real tide between 
stations that it seems preferable to estimate and compare the 
tidal components directly. Moreover, since the low-frequency 
signals are best derived from the tidal residuals, we begin with 
an analysis of the tidal signals. 

Comparison to Tide-Gauge Data: Tidal Regime 

A regional picture of tides covering Lower Cook Inlet and 
the adjoining Gulf of Alaska is shown in Fig. 4 in terms 
of cotidal charts for the largest semidurnal (M 2 ) and diurnal 
(Ki) constituents. These charts are based on data that we have 
extracted from a North Pacific finite-element numerical tide 
model [13]. The semidiurnal M 2 is seen to gain in amplitude 
after entering Cook Inlet, exceeding 200 cm at the mouth 
of Kachemak Bay. The diurnal K\ also grows in amplitude, 
although more weakly. While the charts of Fig. 4 give a useful 
regional view, the resolution of the finite-element model in 
Kachemak Bay is probably too coarse (approximately 4 km) 
to give reliable results in the bay itself, particularly in the upper 
bay. The model does suggest that K\ is fairly constant within 
Kachemak Bay and that Seldovia and PBAY diurnal tides might 
be comparable. It also suggests, however, that the semidiurnal 
tides may differ by several centimeters or more. 



Fig. 4. Cotidal charts of the (left) M2 and (right) Kl tidal coefficients from 
[13]. Amplitudes (in centimeters) are in blue. Greenwich phase lags are in red. 
Black and red circles denote locations of Seldovia and the PBAY GPS site, 
respectively. 

This regional picture of Kachemak Bay tides is gener- 
ally confirmed by historical tide data within the bay proper. 
Space precludes a detailed comparison, but they confirm a 
pattern of nearly constant diurnal amplitudes and phases within 
Kachemak Bay. The semidiurnal M 2 amplitudes increase by 
7 cm from Seldovia to the head of the bay, as do S 2 amplitudes. 
This suggests that the tides at PBAY and Seldovia may differ 
by several centimeters in both M 2 and S 2 , with PBAY being 
larger and only slightly lagged, while the diurnal tides may be 
expected to be very close. 

Based on the Seldovia tide-gauge spectrum, we estimated a 
total of 91 constituents by least squares harmonic analysis. The 
analysis included standard adjustments for nodal and perigee 
modulations of lunar constituents [14]. The variance of the 
residual time series after tidal analysis at Seldovia is 198 cm 2 . 
At PBAY, the residual variance is 386 cm 2 , indicative of higher 
noise in the individual GPS sea-level estimates. For the GPS 
tidal estimates, the standard errors were computed from the 
standard deviations of 12 monthly subset solutions, scaling by 
12 -1 / 2 . (Since Pi and K 2 are unresolved in monthly analyses, 
their standard errors are simply set to those of Oi and S 2 , 
respectively.) These errors are understandably much larger than 
those for the long Seldovia time series but are much less than 
1 cm for even the largest constituents. 

Table I presents a summary compilation of the estimated 
tides. In it are included all diurnal and semidiurnal constituents 
of amplitude 7 cm or greater and all higher order compound 
tides of amplitude 1 cm or greater. Standard errors are also 
included. 

Long-period tidal constituents were also computed but are 
purposely not included in the table to prevent readers from 
thinking they might be legitimate estimates. Only a long mul- 
tiyear time series can provide reliable estimates, and even 
then, they typically have large uncertainties. An analysis of 
one year of Seldovia gauge data, corresponding to the same 
time as the GPS series, yields reasonably good agreement at 
the long-period tidal frequencies, the largest being the annual 
with (H, G ) = (6.8 cm, 227°) for the GPS and (7.4 cm, 233°) 
at Seldovia. The long-term mean of the annual constituent at 
Seldovia is (12.7 cm, 256°). 

A few aspects of Table I are immediately noteworthy. The 
three largest semidiurnal tides are slightly larger in the GPS 
than the Seldovia gauge, and with greater phase lags, and 
this surely owes to real ocean differences between the sites 
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TABLE I 

Amplitude H, Phase Lag G, and Standard 
Error of Selected Tidal Coefficients 



Seldovia gauge 

Peterson Bay GPS 

Tide 

H 

G 

0 

H 

G 

0 


(cm) 

(deg) 

(cm) 

(cm) 

(deg) 

(cm) 

Oi 

34.33 

263.71 

0.04 

33.54 

264.2 

0.42 

Pi 

17.62 

277.90 

0.04 

17.26 

278.7 

0.42 

Ki 

56.10 

279.76 

0.05 

54.63 

278.7 

0.53 

M2 

8.90 

276.81 

0.07 

8.57 

272.8 

0.47 

n 2 

46.37 

298.16 

0.15 

47.31 

299.0 

0.47 

V2 

7.62 

301.77 

0.10 

7.60 

301.2 

0.47 

m 2 

222.27 

324.04 

0.16 

225.43 

324.7 

0.71 

s 2 

82.90 

358.25 

0.08 

83.88 

359.1 

0.72 

k 2 

22.63 

350.89 

0.07 

22.34 

352.5 

0.72 

mo 3 

2.78 

183.21 

0.04 

2.88 

192.8 

0.36 

so 3 

1.37 

209.07 

0.03 

1.31 

211.1 

0.41 

mk 3 

1.88 

201.86 

0.03 

2.10 

202.9 

0.41 

mn 4 

1.09 

40.47 

0.02 

1.48 

45.4 

0.28 

m 4 

3.04 

69.16 

0.03 

3.93 

76.5 

0.28 

ms 4 

2.56 

116.00 

0.02 

3.51 

120.0 

0.48 

m 6 

1.56 

195.77 

0.02 

1.93 

201.1 

0.16 

2MS 6 

1.45 

240.21 

0.03 

1.79 

249.7 

0.28 


as suggested by Fig. 4. The smaller semidiumals agree to 
within their standard errors, except for (the partly nonlinear) /i 2 , 
which displays a slightly larger difference. The largest vector 
difference is 4.1 cm, for M 2 , but this is only 1.8% of the large 
M 2 amplitude. S 2 differs similarly by 1.9%. The difference of 
4.1 cm in M 2 is consistent with what might be inferred from 
Fig. 4 and the historical gauge data. 

A noteworthy feature of the much smaller compound tides 
is that every one shows a larger phase lag at PBAY than at 
Seldovia. Every one, save SO 3 , also shows larger amplitude 
at PBAY. We have no knowledge concerning the generation 
regions for compound tides in the Lower Cook Inlet, but the 
consistent amplitude and phase differences suggest true tide dif- 
ferences at our two locations, with possible wave propagation 
into the bay from the west. 

The case of K\ is of special interest because its period is 
exactly 1 sidereal day, so it is very nearly locked to the GPS 
satellite sampling period. Thus, the K\ frequency can become 
corrupted by any number of errors related to GPS sampling, 
particularly orbit errors [15]. K 2 , locked to half the GPS 
sampling period, can be similarly corrupted. Our results for K\ 
are far from poor and, in fact, are reasonably good with a vector 
difference of 1.8 cm, but they do appear slightly anomalous. 
The difference of 1.8 cm is much larger than corresponding 
differences for the other diurnals 0\ and P\. Moreover, the 
estimated amplitude of K\, 54.6 d= 0.5 cm, is smaller than the 
amplitudes at all the historical stations (not shown), which are 
between 56.1 and 57.2 cm. Our GPS estimate may thus be in 
error by 1-2 cm. 

Not included in Table I is Si, of period of exactly 24.0 h. 
Si is generally a very small tide [16] but of special interest in 
testing new measurement systems because it is often associ- 
ated with thermal errors in instrumentation owing to the daily 
heating and cooling cycle. To compare Si solutions, we use 
our estimates for Seldovia that correspond to the same one-year 
interval as the GPS data. For the tide-gauge data, we obtain 
(77, G) = (1.41 cm, 72.9°); for the GPS, we obtain (1.65 cm, 



Fig. 5. (Top) Daily mean sea levels from the Seldovia tide gauge and PBAY. 
(Bottom) Differences between the top-panel curves. 

74.9°), with a vector difference of only 0.24 cm. Since part of 
this difference could again stem from real ocean differences at 
the sites, we conclude that any errors associated with daily ther- 
mal or radiational effects must be very small in both systems. 

In summary, this analysis of the tidal signals suggests that 
the GPS-based tide estimates appear very reliable and that most 
of the discrepancies with the tide gauge at Seldovia can be 
attributed to real ocean-tide differences between the two sites. 

Comparison to Tide-Gauge Data: Low-Frequency Regime 

In this section, our goal is to quantify the level of daily 
differences and to determine whether any long-term trends are 
evident in the two time series. The Seldovia gauge data used 
in this section are daily means computed using procedures 
identical to those used by the University of Hawaii Sea Level 
Center. This is done by removing a predicted tide from hourly 
data, applying an antialiasing low-pass filter, and computing 
means for each 24-h segment. 

We followed a similar procedure with the GPS time series, 
except for an initial step needed to first form a time series 
of hourly means. This was done from the original unequally 
sampled sea-level estimates by removing a predicted tidal 
signal based on coefficients as in Table I and then applying 
a Savitzky-Golay filter of full width 2.4 h. These filtered 
tidal residuals were then binned into hourly means, and sub- 
sequently, the same procedures as used for the gauge data were 
applied. The one nonstandard aspect of our processing is that 
we have removed from each data set the long-period, as well as 
the short-period, tides; this was done because of small observed 
differences in the long-period tides at the two locations. 

Fig. 5 shows daily mean sea levels for both stations, as well 
as their differences; the rms difference is 2.3 cm. There is no 
hint of any long-period offsets or trends in the differences. As 
before, it is possible that some part of this 2.3 cm owes to true 
ocean differences at the two sites, for example, from differ- 
ences in wind setup. Although we removed long-period tides 
separately from each time series, including annual and semi- 
annual components, which will act to reduce site differences, 
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Fig. 6. Supplemental images for Fig. 1 (left), Fig. 2 (middle), and Fig. 3 (right). See also http://ieeexplore.ieee.org. 


small differences at nontidal periods could remain. We can 
conclude that the GPS technique is measuring daily mean sea 
level at PBAY with at least an accuracy of 2.3 cm and possibly 
somewhat better. 

IV. Conclusion 

We have demonstrated that a geodetic-quality GPS unit sited 
close to the ocean can accurately measure sea level. The benefit 
of a GPS tide gauge is twofold. First, it can measure sea 
level in a global terrestrial reference frame and would thus be 
unaffected by offsets caused by earthquakes and local land sub- 
sidence. A second benefit is that the data from such a GPS unit 
can be used for other purposes, i.e., high-precision scientific 
geodesy/geophysics experiments and as a base station for sur- 
veyors. Although GPS units certainly do fail, most GPS net- 
works operating today recover 99% of the expected data. This 
suggests that GPS tide gauges might be particularly useful in 
low-temperature areas such as Alaska and Antarctica. A true 
side-by-side comparison of a GPS tide gauge and a traditional 
tide gauge is needed to further assess the accuracy of the 
GPS system. 

Appendix 

Fig. 6 provides supplements for Figs. 1-3. 
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